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Some  Comments  on  the  Endochronic  Theory  of  Plasticity 

by 

R.S.  Rivlin 

Center  for  the  Application  of  Mathematics 
Lehigh  University,  Bethlehem,  Pa. 

1 .  Introduction 

In  two  papers  [1,2]  published  in  1971,  Valanis  claimed  to  der¬ 
ive  from  thermodynamic  considerations  three-dimensional  constitutive 
equations  on  the  basis  of  which  the  behavior  of  plastic,  or  visco¬ 
plastic  materials,  subjected  to  small  or  large  deformations  could  be 
modeled.  His  theory  was  most  fully  developed  in  the  case  when  the 
material  is  rate- independent  and  the  deformations  are  sufficiently 
small  so  that  the  classical  infinitesimal  strain  tensor  can  be  used. 
Accordingly,  in  the  present  paper  the  theory  will  be  discussed  in 
this  context. 

In  [1]  and  [2],  Valanis  claimed  for  this  theory  the  merit,  over 
most  other  phenomenological  theories  which  purport  to  describe  the 
elastic-plastic  behavior  of  metals,  that  it  avoids  the  introduction 
of  a  yield  surface.  The  stress  is  expressed  as  a  functional,  of 
specific  form, of  the  history  of  the  strain  -  a  single  such  expression 
being  valid  both  under  conditions  of  loading  and  unloading.  It  also 
expresses  the  stress  explicitly  as  a  functional  of  the  total  strain 
and  avoids  the  separation  of  strain  into  elastic  and  plastic  parts. 

Since  the  thermodynamic  content  of  Valanis' s  argument  is  ques¬ 
tionable,  and  is,  in  any  case,  replete  with  ad  hoc  assumptions,  dis¬ 
cussion  of  it  is  relegated  to  Appendix  A.  There  the  argument  is 
presented  in  a  manner  which  retains  its  essential  features,  but  is 
stripped  of  irrelevancies .  However,  even  when  presented  in  this 
manner,  it  is  evident  that  certain  assumptions  are  unacceptable 
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and  others  are  of  a  quite  arbitrary  character  unsupported  by  physi¬ 
cal  justification.  Accordingly,  in  the  main  body  of  the  paper 
we  discuss  the  constitutive  equation  without  relation  to  this  back¬ 
ground,  purely  from  the  point  of  view  of  its  suitability  as  a  model 
for  such  materials  as  metals. 

The  essential  feature  of  the  constitutive  equation  of  Valanis 
is  that  it  assumes  the  stress  following  some  deformation  history  to 
be  a  linear  isotropic  functional  of  the  strain  increment  history, 
the  kernel  in  this  functional  being  a  scalar  function  of  an 
intrinsic  time  defined  by  an  isotropic  scalar  functional  of  the 
strain  increment  history.  Valanis  uses  the  term  endochronic  theory 
to  describe  a  theory  of  this  type.  Valanis  assumes  a  quite  ex¬ 
plicit  form  for  the  functional  defining  the  intrinsic  time.  He 
justifies  it  by  an  alleged  reasonableness,  which  is  subjective,  and 
by  alleged  agreement  between  the  predictions  of  the  theory  and 
experiment.  It  is  seen  in  §  §  2  and  6  that  even  if  we  take  this 
agreement  at  its  face  value,  it  could  not  establish  the  validity  of 
his  assumed  expression  for  the  intrinsic  time,  since  a  wide  variety 
of  other  expressions  would  provide  equally  good  agreement. 

In  this  paper  we  discuss  the  theory  in  a  slightly  broader 
context,  in  which  the  particular  form  for  the  intrinsic  time  used 
by  Valanis  is  replaced  by  a  somewhat  more  general  one  (see  equation 
5.4)  below).  It  will  be  seen  that  most  of  the  criticisms  which 
can  be  levelled  against  the  theory  of  Valanis  also  apply  to  theories 
based  on  this  more  general  definition  of  intrinsic  time. 

For  rate- independent  materials,  the  strain-history  can 
necessarily  be  parametrized  in  terms  of  the  path-length  covered  by 
the  strain  history  in  9-dimensional  strain  space.  The  intrinsic 
time,  whether  defined  in  the  manner  adopted  by  Valanis,  or  in  the 


more  general  manner  of  the  present  paper,  can  then  be  parametrized 
in  terms  of  this  strain  path  length.  As  a  matter  of  convenience  we 
adopt  this  course  in  the  present  paper,  although  this  is  not 
essential  to  the  arguments  presented. 

In  the  case  of  one -dimensional  deformations,  discussed  in 

§§2-4,  this  strain  path  length  is  taken  as  the  intrinsic  time.  For 

simple  shearing  deformations  this  involves  no  loss  of  generality 

since  the  intrinsic  time  is  necessarily  some  constant  multiple  of 

the  strain  path  length.  This  is  also  the  case  for  simple  extension 

provided  that  we  make  the  assumption,  which  is  made  by  Valanis, 

* 

that  Poisson's  ratio  (defined  as  the  ratio  decremental  transverse 
strain/ incremental  longitudinal  strain)  is  constant. 

In  §2  we  establish  certain  relations  which  could,  at  any  rate 
in  principle,  be  used  to  test  the  applicability  of  the  class  of 
constitutive  equations  considered  to  a  particular  material  for  one¬ 
dimensional  deformations. 

In  §3  we  discuss  the  relation  between  the  incremental  moduli 
for  loading  and  unloading  from  a  state  in  the  plastic  regime.  It 
has  been  pointed  out  that  the  constitutive  equation  adopted 
by  Valanis  yields,  in  disagreement  with  experimental  results  for 
metals,  a  value  for  the  unloading  modulus  which  is  very  different 
from  the  initial  modulus  for  infinitesimal  strains  from  the  un¬ 
deformed  state.  We  find  that  this  is  not  necessarily  the  case  in 
the  context  of  the  more  general  theory  discussed  in  the  present 
paper.  In  §4,  we  discuss  the  dissipation  in  infinitesimal 
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This  assumption,  which  was  criticized  by  Lee  [3]  is  discussed 
in  §  7 . 


unloading-loading  and  loading-unloading  strain  cycles  starting  at 
an  arbitrary  strain  which  has  been  reached  by  a  monotonic  loading. 

We  find,  as  did  Sandler  [4]  in  a  more  restricted  context,  that  the 
dissipation  must  necessarily  be  negative  in  either  one  or  the 
other  of  these  cycles  and  accordingly,  the  material  will  be  unstable. 

The  three-dimensional  theory  is  introduced  in  §  5 .  It  is 
shown  in  §6  that  for  deformations  consisting  of  successive  simple 
extension  and  simple  shears,  carried  out  in  discrete  time  inter¬ 
vals,  any  definition  of  intrinsic  time  in  the  broader  class  envisaged 
in  the  present  paper  can  lead  to  the  same  prediction  for  the  stress 
as  that  used  by  Valanis,  if  appropriate  values  are  given  to  the 
adjustable  constants.  Accordingly,  agreement  between  theory  and 
experiment  for  deformations  of  this  type  cannot  be  used  to  estab¬ 
lish  the  predictive  value  of  Valanis 's  constitutive  equation  for 
more  general  deformations.  Attention  is  drawn  in  §8  to  a  further 
general  criticism  of  constitutive  equations  based  on  the  endo- 
chronic  concept.  It  is  shown  that  if  a  strain  is  reached  by  two 
paths  which  differ  only  very  slightly,  the  associated  stresses  may 
be  very  different,  i.e.  the  relation  between  stress  and  strain 
history  is  not  continuous  in  the  sense  of  the  supremum  norm. 

In  [5]  Valanis  presented  a  revised  theory  in  order  to  meet 
the  criticisms  which  had  been  levelled  at  the  earlier  theory. 

This  is  discussed  briefly  in  §9.  He  again  motivates  the  theory 
by  thermodynamic  arguments  which  are  analogous  to,  and  open  to  the 
same  criticisms  as,  those  used  to  motivate  the  earlier  theory. 
Accordingly,  we  do  not  discuss  them  explicitly  in  the  present 
paper.  While  this  second  theory  is,  in  fact,  not  open  to  most  of 
the  criticisms  to  which  the  earlier  theory  is  subject,  it  does  not 
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possess  the  two  major  and  somewhat  revolutionary  advantages  which 
were  claimed  for  the  latter  -  that  it  describes  the  elastic-plastic 
behavior  of  metals  without  the  introduction  of  a  yield  surface  and 
without  the  separation  of  the  strain  into  elastic  and  plastic 
parts.  It  differs  from  previous  theories  of  plasticity  only  in  the 
manner  in  which  strain-hardening  affects  the  yield  surface  and  the 
relation  between  stress  and  incremental  plastic  strain.  No  evidence 
is  adduced  to  demonstrate  that  it  provides  a  more  accurate  des¬ 
cription  of  the  actual  behavior  of  metals  than  do  other  constitu¬ 
tive  equations  that  have  been  suggested. 


2 .  One-dimensional  theory 

We  consider  small  uniaxial  deformations  of  a  rate-independent 
material.  Let  e  denote  the  strain.  Let  Z  be  the  length  of  the 
strain-path,  thus: 


d Z  =  | de  j  ,  (2.1) 

with  Z  =  0  in  the  undeformed  state. 

We  now  consider  a  deformation  in  which  t  =  0  initially  and  in¬ 
creases  monotonically  to  e  ,  then  decreases  monotonically  to 
c2  ,  increases  monotonically  to  e 3  ,  and  so  on,  the  final  rever¬ 
sal  taking  place  at  e  =  .  We  denote  the  final  value  of  t  at 

which  the  stress  is  measured  by  E  .  We  note  that  if  y  is  odd  (even 
the  last  reversal  will  be  followed  by  a  decrease  (increase)  in  e  . 

Let  be  the  value  of  Z  corresponding  to  the  ath  rever¬ 

sal  of  the  strain,  and  let  L  be  the  value  of  l  corresponding  to 
the  final  value  E  .  It  is  easily  seen  from  (2.1)  that 

Ct-1 

1=2  l  (-1)6"1  efl  +  (-l)a~1  e  ,  (2.2) 

a  g“x  S  v  a  ’  v  J 

and 

l  •  2  l  C-l)8"1  ea  *(-1):j  E  .  (2.J) 

B=1  3 


We  now  regard  e  as  a  function  of  Z  and  make  the  assump¬ 
tion  that  the  stress  a  is  given  by 


a 


L 

f(L,*)de(£)  , 

0 


(2.4) 


where  f  is  a  positive  function  of  the  indicated  arguments. 
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Noting  that  for  L  .<  i  <  l  , 

a-l  a 


de(£)  =  d£  (a  odd)  and  de(£)  =  -d l  (a  even)  ,  (2.5) 


we  obtain  from  (2.4) 


o-2  £  C-D01'1  g(L,£  )  +  (-1)M  g(L,L)  , 

a=l  a 


where 


(2.6) 


g(L,f)  = 


f(L,f)d l  . 


(2.7) 


We  can  rewrite  (2.6)  as 


a-l 


(-l)a-1  a(L,£a)  +  o(L,L)  (u  even) 


a  = 


l  (-l)a~1  a(L ,Z  ) 
a-l  a 


(u  odd)  , 


(2.8) 


where 


a(L,£a)  =  2g(L,£a)  -g  (L,L)  , 
a(L,L)  -  g(L,L) 


(2.9) 


We  note  that  a(L,£a)  is  the  stress,  when  the  length  of  the  strain 
path  is  L  ,  in  a  deformation  in  which  the  strain  increases  mono- 
tonically  from  zero  to  i and  then  decreases  monotonically  until 
the  total  length  of  the  strain  path  is  L  .  Also,  o(L,L)  is  the 
stress,  at  strain  path  length  L  ,  in  a  deformation  in  which  the 
strain  increases  monotonically  from  0  to  L  . 

From  (2.8)  it  is  easily  seen  that 
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a(L,*lf...,£u)  -a(L,£1>...>£u_1)  =  (-l)M[a(L,L)  -o(L^u)],  (2.10) 

where  a(L,£1, . . . ,1  )  denotes  the  stress,  at  strain  path  length 
L  ,  resulting  from  a  deformation  with  u  reversals  of  strain  at 
path  lengths  . 

From  (2.10)  or  (2.8),  we  easily  obtain  the  relations 


ucl,^ . ia)  -ocl,^ . £0_2)  - 


Conversely,  by  giving  a  the  values  3,5,...,y  (odd)  in  (2.11) 
and  adding  the  resulting  equations,  we  can  recover  (2.8)2  and,  by 
giving  at  the  values  2,4,...,y  (even)  and  adding  the  resulting 
equations,  we  can  recover  (2.8)^  .  In  deriving  the  second  of  these 
results,  we  note  that  when  a  =  0  ,  a (L ,1^) =a (L ,L)  . 

As  an  example  of  the  implications  of  equations  (2.10)  and  (2.11), 
we  consider  a  deformation  in  which  the  strain  increases  from  zero 
to  e  and  then  undergoes  n  cycles  of  deformation  between  strains 
e  and  e2  (s^s^)  •  In  this  case 


ia  =  ex  +  (a-l)(e1-s2)  ,  L  =  e1+  2n(el-e2)  . 


(2.12) 


Taking  ct  =  2n  +  1  in  (2.11),  we  obtain 


,  t  ^ ,  . . . ,  -^2n+l^  ^  ’  ^1  ’  *  ’ '  ’  ^2n  - 1^  ”  cr^’^2n+l^  'a^’^2n^  *  (2.13) 

We  note  that  the  a’s  in  (2.13)  have  the  following  interpretations: 
o(L,£1,  . . .  ,£2n+i)  is  the  stress  at  strain  z-^  resulting  from  an 
initial  deformation  to  strain  followed  by  n  cycles  between 

strains  e  and  z2  (e2<e1)  5  a(L ,^1 , . . . »^2n-l^  is  the  stress  at 


strain  e1+  2(e1-e2)  resulting  from  an  initial  deformation  to 
strain  e1  followed  by  n-1  cycles  between  strains  e1  and  e2 
and  then  a  further  increase  in  strain  to  2(e1-e2)  ; 

a(L,£7n+1)  =  a(L,L)  is  the  stress  at  strain  e^+  2n(e1-e2)  resul¬ 
ting  from  a  monotonic  deformation  to  this  strain; 

a(L,f2n)  is  stress  at  strain  e1+  2  (n-1)  (e1-e2)  resulting  from 

a  monotonic  deformation  to  strain  e1+  (2n-l)(e1-e  )  followed  by  a 
decrease  in  the  strain  to  e1+  2(n-l)(e  -e  )  . 

Suppose  now  that  we  do  not  assume  that  the  stress  is  given  by 
an  expression  of  the  form  (2.4).  As  before,  we  consider  that  ini¬ 
tially  the  strain  increases  and  this  increase  is  followed  by  y  rever¬ 
sals  of  the  strain  at  strains  e  , e  ,  ...,e  .  Then,  at  l  -  L  ,  the 

1  2  y 

stress  must  be  a  function  of  L  and  of  e  ,e  ,...,e  ,  and,  hence 

12  y 

from  (2.5),  of  L ,  l±,  .  We  write 

a  =  a(L,fl,  .  .  .  ,£y)  .  (2.14) 

We  note  that  if  y  =  1,  the  relation  (2.10)  is  satisfied  identi¬ 
cally.  (We  bear  in  mind  that  if  y  =  1,  o(L ,£,..., £^)=a(L ,L) . ) 

It  is  easily  seen  that  if  the  function  cr(L , =  a(L,f1)  ,  say, 

is  specified,  and  we  take 

£(L,£)  =  7  Jl  toCL,^]  (2.15) 

in  (2.4),  we  obtain 

ff( L,f1)  =  \  [2a(L,f1)  -  a (L , 0)  -  a(L,L)]  .  (2.16) 

a(L,0)  gives  the  specified  dependence  of  stress  on  path  length 
for  monotonic  decrease  of  strain  from  zero  to  -E  .  Provided  that 
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a(L,0)  =  -a(L,L)  , 


(2.17) 


we  obtain  from  (2.16) 


a(L,£1)  =  a(L,£1)  .  (2.18) 

The  condition  (2.17)  is  the  condition  that  the  prescribed  stress 
merely  changes  signs  if  the  prescribed  strain  history  is  replaced 
by  its  negative.  That  such  a  condition  is  necessary  if  the  pre¬ 
scribed  stress  is  to  be  expressible  in  the  form  (2.4)  is  evident. 

We  thus  conclude  that  if  the  stress  at  l  =  L,  following  a  single 
reversal,  changes  sign  when  the  strain  history  is  replaced  by  its 
negative,  then  this  stress  can  be  expressed  in  the  form  (2.4). 

We  shall  now  show  that  if,  further,  the  specified  stress  at 
i  -  L  following  y  reversals  of  strain  is  an  odd  function  of  the 
strain  history  and  satisfies  the  relations  (2.10)  for  u  =  l,2,...,y, 

then  the  specified  stress  can  be  expressed  in  the  form  (2.4)  with 

f(L,£)  given  by  (2.15). 

We  prove  this  result  by  induction.  Suppose  that  <j( L, 

Za  )  with  a  odd  is  given  by  (2.4),  i.e.  by  (2.6)  with  y  =  a-1, 
thus 

a  =  2[g(L,£1)  -  g(L ,l2)  *  ...  -g(L,fa_l)]  +  g(L,L).  (2.19) 

We  assume  that  (2.10)  with  y  =  a  is  valid.  Then,  we  obtain, 
from  (2.19),  (2.9)  ,  and  (2.10)  with  y  =  a  , 

a(L,£1,...,fa)  =  2[g(L,f1)-g(L,f2)  +  ...  +  g(L,£a)]  -g  (L ,  L) .  (2  . 20) 
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It  follows  that  o(Lf£  , . . . ,1  )  is  given  by  (2.4)  with  f(L,£) 
given  by  (2.15).  An  analogous  argument  is  applicable  to  the  case 
when  a  is  even.  Since  we  have  already  seen  that  the  proposed 
theorem  is  valid  when  a  =  1  ,  it  is  valid  for  all  a  . 

It  is  suggested  that  an  experimental  test  of  the  validity  or 
non-validity  of  the  constitutive  equation  (2.4)  might  well  be  based 
on  the  relations  (2.10). 


> 


-  12  - 


.  The  incremental  modulus  -  uniaxial  deformation 


In  this  section,  we  again  consider  a  uniaxial  deformation  in 
which  the  stress  is  given  in  terms  of  the  strain  history,  regarded 
as  a  function  of  strain  path  length  Z,  by 


rL 

a(L)  =  £(L,Z)de(l)  , 


(3.1) 


where  L  is  the  value  of  Z  at  which  the  stress  is  measured. 

We  consider  a  strain  history  which  consists  either  of  a  mono- 
tonically  increasing  strain,  or  of  a  sequence  of  increasing  and  de¬ 
creasing  strains.  We  assume,  however,  that  the  strain  E  corres¬ 
ponding  to  path  length  L  is  positive  and  is  reached  finally 
through  increasing  strain.  We  calculate  the  incremental  moduli  at 
l  =  L  when  the  strain  is  increased  by  an  infinitesimal  amount  AE . 
Let  Acr  be  the  corresponding  increase  in  stress.  If  AL  is  the 
increase  in  Z  corresponding  to  the  strain  increase  AE  ,  we  have 


AL  =  I  AE 


(3.2) 


We  denote  by  u+  and  u  the  incremental  moduli  correspond¬ 
ing  to  AE  positive  and  negative  respectively.  Then,  from  (3.1), 


La  _  Act 
l+  ~  AE  AL 


f  (L  ,  L)  +  jj;  f[L,Z)dZ 


(3.3) 


u-  ■  If  ■  'lr  '  f(L-L)  "f 


We  note  that  when  L  =  0  , 


f 
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yQ  is  the  incremental  modulus  at  zero  deformation.  From  (3.3) 

y+  +  y_  =  2f(L,L)  .  (3.5) 

Valanis  [1,2]  made  the  particular  choice  of  f(L,£): 


f(L,£) 


Eo  *  EiirrsL 


(3.6) 


where  EQ,  E  ,  3  and  n  are  positive  constants.  In  this  case, 


y0  =  f (L,L)  =  Eq  +  E1  , 
and  (3.5)  becomes 


(3.7) 


y+  +  u_  =  2yQ  . 


(3.8) 


For  metals  it  is  found  experimentally  that  u_  is  approxi¬ 
mately  yQ  and  since,  in  the  plastic  regime,  p+  is  very  much 
less  than  yQ  ,  the  relation  (3.8)  is  not  valid.  This  was 
adduced  as  an  argument  against  the  endochronic  theory  ad¬ 
vanced  by  Valanis  in  [1,2]  and  led  the  latter  to  modify  his  theory. 
It  is  evident  that  the  same  criticism  would  apply  to  any  theory  in 
which  f(L,£)  has  the  form 


f(L ,£)  =  C1 


^TLj  +  c2 


<t>2(£) 


r 

+  Cv  TJXT  9 


(3.9) 


where  C1,...,Cv  are  constants,  so  that  f(L,L)  =  f(0,0)  .  How¬ 
ever,  it  may  not  apply  for  other  forms  of  f(L,£).  Also,  this 
criticism  does  not,  of  course,  preclude  the  applicability  of  one¬ 
dimensional  constitutive  equations  of  the  type  considered  to  other 
materials  than  metals. 
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4 .  Dissipation  in  cyclic  deformation 
(a)  Unloading -loading 

We  now  consider  the  material  to  undergo  a  monotonically 
increasing  uniaxial  strain  to  strain  e1  .  Thereafter,  the 
strain  is  decreased  to  e2  and  then  increased  to  e1  .  We 
shall  calculate  the  dissipation  V  ,  per  unit  volume,  in  this 
cycle  of  deformation. 

Let  £  ,  £2  and  L  denote  the  strain  path  lengths  at 
the  first  and  second  strain  reversals  and  at  the  final  strain 
e1  respectively.  Then, 

^1  =  £1  ’  l2  =  2ei  '  £2  *  L  =  3ei  ‘  2e2  *  (4-l) 

Let  l  be  the  value  of  the  strain  path  length  at  a  generic 
point  on  the  cycle  at  which  the  strain  is  e  .  Then, 


d l  =  de  (strain-increasing)  and  d l  =  -de  (strain-decreasing) . 

(4.2) 

From  (2.4)  and  (4.2)  the  stress  at  a  generic  point  of  the 
cycle,  corresponding  to  path  length  l  ,  is  given  by 


ff(i,y  =  2g(£,£1)  -g  (£,£) 

=  2g(£,£1)  -2g(£,£2)  +  g (£,£), 

where  g (£,C)  is  defined  by  (cf.(2.7)) 

rC 


g(£»C) 


f (£,5)d5  . 


0 


strain -deer easing 

(4.3) 

strain-increasing  , 


(4.4) 


The  dissipation  V  ,  per  unit  volume,  in  the  cycle  is  given  by 
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V 


aUU^dZ  + 


o(l,Z1,l2')dl  . 


(4.5) 


We  now  make  the  assumption  that  the  amplitude  of  the 
deformation  is  sufficiently  small  so  that  terms  of  second  deg¬ 
ree  in  (£-£  )  may  be  neglected  in  comparison  with  those  of 
first  degree.  Then,  we  may  write,  from  (4.4), 


g U,l±)  =  g U1,l1)  +  C^-^g!  , 

=  g(^1,£1)  +  (l-Z1)g1  +  U2-l1)g2  ,  (4.6) 

g(£,£)  =  g U1,l1)  +  (£"e1)g1+  U-l±)g2  , 

where  gx  and  g2  are  defined  as  the  values  of 

^  and  9*d£'E’^  respectively  when  l  =  l1  and  £  *  l1  . 
Introducing  (4.6)  into  (4.3),  we  obtain 
0(1,1,)  =  g (t  ,1  )  +  U-l1)g1  -  U-£1)g2  , 

(4.7) 

0(l,llfl2)  =  g^]_,f1)  +  U'l1)g1  +  (■^'^■Z.tg)  gg  . 

We  now  introduce  (4.7)  into  (4.5)  and  carry  out  the  integra¬ 
tions  to  obtain,  with  (4.1), 

V  s  glCV£2)2  *  (4.8) 

For  the  particular  expression  for  f(£,5)  used  by 
Valanis  (cf.(3.6)),  viz., 

•  Eo  ‘  Ei(  T^H)""1 


(4.9) 


where  EQ ,  E^,  8  and  n  are  positive  constants,  we  have 
from  (4.4)  and  the  definition  of  g  : 


*1 


E^n-l) 

n 


1 

(i^e^1)n 


(4.10) 


Loading -unloading 

We  now  consider  the  material  to  undergo  a  monotonically 
increasing  uniaxial  strain  to  e  .  It  then  undergoes  a  cycle 
of  strain  in  which  the  strain  is  increased  to  e2  >  e1  and 
then  decreased  to  e1  again.  The  strain  path  lengths 
Z^,  Z2  and  L  are  now  given  by 

l±  =  ,  l2  =  e2  ,  L  *  2e2-e1  .  (4.11) 

If,  as  before,  l  is  the  value  of  the  strain  path 
length  at  a  generic  point  on  the  cycle  at  which  the  strain  is 
e  ,  we  have,  as  before,  the  relations  (4.2).  The  stress  is 
given  by 


ail, l )  =  g (1,1)  strain- increasing 

a{l,Z2)  =  2g (Z,Z2)  -g (Z,Z)  strain-decreasing  . 


(4.12) 


The  dissipation  V  in  the  cycle  is  now  given  by 


o(Z,Z)dZ  - 


aU,l2)dZ 


(4.13) 


Again,  with  the  small  amplitude  assumption,  we  obtain 
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gU,D  -  *  U-LJSz  » 

S(*.*2)  =  g(i1,£1)  +  (^-^1)g1  +  (V£1)s2  ' 


(4.14) 


From  (4.12),  (4.13)  and  (4.14)  we  obtain,  with  (4.11), 


P  s  -g1(e2‘ei^2  • 


(4.15) 


If  f(£,C)  is  given  by  (4.9),  we  find,  as  before,  that  gx 
is  given  by  (4 . 10) . 

(c)  Discussion 

We  note  from  (4.4)  that 


gl 


3gU,q 

it 


I, 


i 

Sl 


d5  . 


(4.16) 


Then,  gx  *  0  if  l^  ■  0  and,  from  (4.8)  and  (4.15),  we  have 
V  *  0  ,  whether  the  infinitesimal  cycle  of  deformation  is  the 
loading-unloading  or  the  unloading-loading  cycle.  If  f  is  a 
monotonically  decreasing  function  of  l  ,  then  gx  is  nega¬ 
tive  and,  for  l^  ft  0  ,  V  <  0  for  the  unloading-loading 
cycle  and  V  >  0  for  the  loading-unloading  cycle.  This  is 
the  situation  which  prevails  in  the  case  when  f 
takes  the  particular  form  (4.9)  proposed  by  Valanis.  These 
results  imply  that  the  materials  considered  have  no  finite 
elastic  range  as  was  pointed  out  by  Lee  [3],  Also,  the  fact 
that  V  <  0  for  the  unloading - loading  cycle  implies,  was 
pointed  out  by  Sandler  [4],  that  the  material  modeled  is 
unstable . 


5. 


Three-dimensional  theory 

We  consider  the  material  to  undergo  sufficiently  small  defor¬ 
mations  so  that  the  deformation  in  an  element  can  be  described  by 
the  history  of  the  infinitesimal  strain  matrix  e(t)  3  llejj  (t)  || 
referred  to  a  rectangular  cartesian  coordinate  system  x  .  Let 
o  3  ||cr^j  ||  denote  the  stress  matrix  at  time  T  ,  referred  to  the 
system  x  . 

The  length  of  the  strain  path  at  time  t  in  9-dimensional 

* 

strain  space  is  given  by 

<U(t)  3  { tr [dc  (t)  ,  £(0)  =  0  , 

i .  e . 

It t)  3  ^  (tr[de(t)  j2}*5  . 

0 

We  introduce  the  notation 

L  3  £( T)  .  (5.3) 

It  has  been  shown  by  Pipkin  and  Rivlin  [6]  that,  for  an 
isotropic  rate-independent  material,  the  stress  corresponding  to 
strain  path  length  L  is  an  isotropic  tensor  functional  of  the 
strain  history  regarded  as  a  function  of  l  .  We  shall  examine 
the  properties  of  a  particular  class  of  such  relations  which  in¬ 
cludes  that  used  by  Valanis  [1,2]  as  a  special  case. 


(5.1) 

(5.2) 


*  The  strain  path  lengths  l  and  L  introduced  in  this  section 
do  not  reduce  in  the  uniaxial  case  to  those  introduced  in  §2. 
They  can,  however,  be  simply  related  to  them. 
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We  define  an  isotropic  scalar  functional  ?  of  e(£)  by 

d?  =  d s(£)  =  4>[deU)]  ,  c(0)  =  0  ,  (5.4) 


where  d>  is  a  positive  definite  isotropic  scalar  function  of 
de(£)  ,  homogeneous  of  degree  unity  in  the  latter.  From  (5.4) 


5  = 


♦[dett')]  . 


(5.5) 


We  call  c  the  intrinsic  time.  We  note  that  since  <p  is  posi¬ 
tive  definite,  q  increases  monotonically  with  L.  Since  $  is  an  iso¬ 
tropic  function  of  de(£)  ,  it  must  be  expressible  as  a  function 
of  the  elements  of  an  isotropic  integrity  basis  for  de(£)  .  This 
may  be  chosen  as  1^,  Ig,  1^  defined  by 

l±  =  tr  [de  (£)  ]  ,  I2  *  tr  [de  (£)  ] 2  ,  I3  -  det[de(£)]  .  (5.6) 


We  adopt  the  notation 


Z  =  C  CL)  . 


(5.7) 


We  now  assume  that  the  stress  a  corresponding  to  strain  path 
length  L  is  given  by  an  expression  of  the  form 


a=2 


u(Zfc)d  e(£)  +  6 


X(Z.O  tr [de (£) ]  , 


(5.8) 


0  0 
where  <$  denotes  the  unit  matrix. 

Valanis  adopts  the  particular  form  for  <J>  : 


4> [ de  (£)  ]  =  (kxI l  *  k2I2)ls  , 


(5.9) 


where  k1  and  k2  are  positive  constants,  which  may  vary  from 
material  to  material.  He  gives  no  reason  for  making  this  particu¬ 
lar  choice  beyond  the  statement  "it  appears  logical  to  define  5 
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6 •  Simultaneous  simple  extension  and  simple  shear 

We  now  suppose  that  the  material  is  subjected  to  a  simple 
extension  in  the  x^-direction  of  a  rectangular  cartesian  coordin¬ 
ate  system  x  and  a  simple  shear  for  which  the  direction  of  shear 
is  the  Xg-direction  and  the  plane  of  shear  is  the  x^x^-plane. 

We  adopt  the  notation 

e  =  eU)  =  e1;LU)  ,  k  =  <(£)  =  e12U) 

and 

E  =  e(L)  ,  K  =  k (L)  , 

where  Z  -  L  at  the  instant  at  which  the  stress  is  measured. 

In  order  to  be  specific  we  shall  suppose  that  the  initial 
extension  is  positive  and  that  successive  reversals  of  the  exten¬ 
sion  occur  at  e  *  e1,e2> • • • >ep  .  We  also  assume  that  the  initial 
shear  is  positive  and  successive  reversals  occur  at  K  =  <1,ic 
It  is  assumed  that  the  shearing  and  extensional  deformations  may 
occur  in  any  order,  but  take  place  in  disjoint  time  intervals. 

Let 


il 

» 

when 

e  =  ea  (a=1” 

• . »u)  , 

1  =  7 

0 

when 

K  -  <g  (8=1,. 

.  . , v)  . 

For  a  simple  extension  in  the  x  -direction  only  o11  is  non¬ 
zero  and,  from  (5-8),  for  the  simple  shearing  deformation  in  the 
12-plane,  only  a12(=a21)  *s  non‘zer°. 

We  assume,  with  Valanis,  that  Poisson's  ratio  for  the  material 
is  a  constant  u  ,  say,  independent  of  the  deformation  to  which 
the  material  is  subjected.  (We  shall  discuss  this  restriction  later 


(6.1) 

(6.2) 
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in  i 7) .  Then, 


tr  e  =  (l-2w) e  . 


(6.4) 


IVe  have,  from  (5.2),  (6.3)  and  the  assumption  that  Poisson's  ratio 
is  constant 


l  3  (1  +  ^2)h\  l  2  (-l)a-1e  +  (-l)ae 


+  2 h\  ^  2  (-1  )ct-1<a  +  (-l)a<l  , 


(6.5) 


where  0  is  the  number  of  reversals  of  extensional  strain  and  y 
is  the  number  of  reversals  of  shear  strain  prior  to  the  strain-path 
length  l  being  attained. 

Since  $  is  homogeneous  of  degree  unity  in  de(£)  and  is 
positive  definite,  for  the  particular  class  of  deformations  con¬ 
sidered  in  this  section 


k  [de(£)|  for  simple  extension 
-  _  (6.6) 
k2Jd<(£)|  for  simple  shear  , 


where  k1  and  k2  are  positive  constants,  which  depend  on  the 
particular  choice  of  <J>  .  We  note  from  (6.5)  that 


| de (£) |  -  (l  +  2u>2) d l 

|d<(£) |  =  d l 


(6.7) 


From  (6.4)  and  (5.6),  we  have 


IL  -  (l-25)de,  I2  *  (l+2u2)  (de) 2 ,  I3 
Ix  -  I3  -  0  ,  I2  -  2 (d<) 2 


^  p  3 

u)  (de)  for  simple  extension 

(6.S) 

for  simple  shear. 
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In  the  particular  case  (5.9)  considered  by  Valanis,  we 
obtain,  with  (6.7)  and  (6.8), 

v>  =  {k1(l-2w)2+  k2(l  +  22“)}2  (1  +  222)-Js  d  l  for  simple  extension 
<p  -  k2  dl  for  simple  shear  . 

Introducing  (6.7)  into  (6.6)  and  comparing  the  resulting  expres¬ 
sions  with  (6.9),  we  obtain 


kx  =  (kx (1-22) 2  +  kp(l+222)}^ 
k2  =  (2k2)^  . 

As  another  illustration,  suppose 
4>  =  [alj  +  bl2l2  +  cI3IL  +  el2]k  . 


(6.10) 


(6.11) 


where  a,b,c,e  are  positive  constants  and  I1,  I2,  I3  are  de¬ 
fined  by  (5.6).  Then,  from  (6.11),  (6.7)  and  (6.8) 


$ 


[a (1 -  22) ^  +  b (1 -  22) 2 (1  +  222) 

+  e (1+222) 2]k  (1+222) 


e^  d l 


+  c22(1-22) 

-h 

d l  for  simple  extension 

(6.12) 

for  simple  shear. 


Comparing  (6.12)  with  the  expressions  for  <J>  obtained  from  (6.6) 
and  (6.7),  we  have 


kx  =  [afl-22)4  +  b  (1  -  22)  2  (1  +  222)  +  c22(l-22)  +  eU+222)2]*5  , 
k2  =■  (4e)k  . 


(6.13) 
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Iv'e  have  seen  that  any  b  ,  which  is  a  function  of  I1,  Ip>  i 
and  is  homogeneous  of  degree  1  in  de/dt  ,  leads  to  the  expres¬ 
sion  (6.6)  for  simple  extensions  and  simple  shears  of  the  types 
considered.  It  follows  that  we  cannot  distinguish  between  differ 
ent  forms  of  <p  by  experiments  in  which  the  deformation  is  a 
sequence  of  such  simple  extensions  and  simple  shears.  In  particu 
lar,  the  experiments  of  Mair  and  Pugh,  which  involve  the  super¬ 
position  of  simple  extensions  on  simple  torsions  of  thin  tubes, 
cannot  be  used,  as  Valanis  has  done,  to  establish  the  validity  of 
the  particular  form  (5.9)  for  cp  .  Agreement  of  his  theory  with 
such  experiments  could,  at  best,  lead  only  to  the  conclusion  that 
the  theory  is  not  inconsistent  with  the  experiments. 


it 
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7 .  Poisson's  ratio 

It  is  assumed  by  Valanis  in  [1,2]  that  Poisson's  ratio  for 
the  materials  he  considers  is  constant.  An  assumption  as  radical 
as  this  should  be  tested  independently  for  any  material  to  which 
the  theory  is  applied.  This  could  be  done,  at  any  rate  in  prin¬ 
ciple,  by  measuring  the  change  in  volume,  or  the  lateral  contrac¬ 
tion,  of  a  rod  of  the  material  when  subjected  to  simple  extension. 

Alternatively,  it  might  be  done  by  making  simultaneous 
measurements  of  tensile  and  shearing  force  when  a  thin  cylinder  of 
the  material  is  simultaneously  subjected  to  simple  extensional  and 
shear  strains  which  are  increased  proportionately.  For  simultane¬ 
ous  monotonically  increasing  simple  extensions  and  simple  shears,  we 
find  from  (5.8)  that,  whether  or  not  Poisson's  ratio  u>  is  constant 


|  [2u(Z,0  +  (l-25)*(Z,C)]deU)  , 

0  rL 

„  =  [-2Su(Z,c)  +  (l-2(3)X(Z,0]de(£) 


2u(Z,?)d<(f)  ,  a23  =  a31  =  0 


where  e(£)  and  <(£)  denote  the  extensional  and  shear  strains 
respectively.  In  (7.1),  l  is  given,  from  (5.1),  by 


(7.1) 


d l  =  [  (l+2w2)  (de)  2  +  2  (die)  2 ] ^  , 


(7.2) 


with  l  =  0  when  the  material  is  undeformed.  Also,  from  (5.6), 
we  have 

IL  =  (1-  2w)  de  (£)  ,  I2  =  (l  +  2ui2)  [de  (£)  ]  2  +  2[d  <Uj)2  , 

I,  -  i[de(f)]{w[de(f)]2  +  [d^(f)]2}  . 


(7.3) 
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d;  is  a  positive  function  of  I  ,  I  I , ,  homogeneous  of  degree 
unity  in  de(f)  and  d<(2),  and  ;  =  0  when  1=0. 

It  follows  from  (".1)  that,  if  2  is  constant, 


axl  =  2(1+2) 


uC,;)d e(£)  , 


,L 

0 


o12  =  2  u(Z,;)d<U)  . 


(7.4) 


Now,  if  the  extensional  and  shear  strains  are  increased  propor¬ 
tionately,  so  that 


=  xe(£)  > 

where  y  is  a  constant,  we  have,  from  (7.2)  and  (7.5), 

de(£)  =  (l+222+  2X2)~h  d l  , 
d<(f)  =  x(1  +  222+  2x2)_}5  df  . 


(7.5) 


(7.6) 


Also,  from  (7.4)  and  (7.6),  we  have 


rxl  =  2(1  +  2)  (1  +  222+  2x2)  ^  j  y(Z,c)df 

0 


a12  =  2(1+22  +  2y  )“'s  x 


u(Z,g)d£  . 


Thus , 


ail/al2  =  (1+35/x  •  (7-3) 

The  assumption  that  Poisson's  ratio  is  constant  has  been  criti¬ 
cized  by  Lee  [3]  insofar  as  the  applicability  of  the  theory  to 


27 


metals  is  concerned.  For  most  metals,  Poisson's  ratio  in  the 
elastic  range  is  about  1/4.  This  implies  that  the  ratio 
bulk  modulus/Young's  modulus  ~  2/3.  It  is  well-known  that  when 
the  plastic  strain  in  a  metal  is  large  compared  with  the  elastic 

A 

strain,  the  bulk  modulus  is  far  greater  than  the  Young's  modulus  . 

It  is,  however,  possible  that  for  other  materials,  particu¬ 
larly  those  which  may,  with  good  approximation,  be  regarded  as  in¬ 
compressible,  the  assumption  of  constant  Poisson's  ratio  will  be 
valid . 


*  Here  we  are  using  the  term  "Young's  modulus"  in  the  sense  of 
(tensile  force  per  unit  area/extensional  strain)  notwithstanding 
that  the  material  is  plastically  deformed. 
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8 .  A  continuity  consideration 

We  now  consider  a  deformation  which  consists  of  alternate 
infinitesimal  simple  extensional  strains  of  amount  Ae  and  in¬ 
finitesimal  simple  shearing  strains  of  amount  Ak  .  As  in  §6,  we 
consider  the  extension  to  be  parallel  to  the  x^axis  of  a  rectan¬ 
gular  cartesian  coordinate  system  x  and  the  shear  to  be  in  the 
x2-direction  and  the  x^Xg-plane.  We  consider  the  tensile  and 
shearing  stresses  after  v  such  simple  extensions  and  simple 
shears.  We  denote  the  resultant  extensional  strain  by  E(=vAe) 
and  the  resultant  shear  strain  by  K(=vA<)  .  The  value  of  t,  when 
the  extensional  strain  is  £  and  the  shearing  strain  is  <  is 
given  by  (cf . (6.6)) 


S  =  kxe  +  k2<  .  (8.1) 

Since 

< / £  =  K/E  ,  (8.2) 

equation  (8.1)  can  be  rewritten  at 

5  -  (*!  +  |  k2)s  -  Cf  +  k2)<  .  (8.3) 


We  shall  assume  that  r  has  the  form  proposed  by  Valanis  and 
given  by  (5.5)  and  (5.9).  Then,  k  and  k2  are  given  by  (6.10). 

With  the  assumption  that  Poisson's  ratio  is  constant,  the 
tensile  and  shear  stresses  are  obtained  from  (7.4)  as 


of11  =  2(l+<2) 


fE 


0 


q(Z,c)de  , 


°12  =  2 


M  (Z  ,  C)  die  , 


(8.4) 


0 


where  £  is  given  by  (8.3)  and 


Z  =  knE  +  kgK  . 


(8.5) 


We  now  suppose  that  instead  of  increasing  the  strains  to  their 
final  values  in  a  stepwise  fashion,  we  increase  them  proportion¬ 
ately.  Then,  using  the  expression  for  £  given  by  (5.5),  (7.3) 
and  (5.9)  we  have,  with  (8.2)  and  (6.10), 


(8.6) 


Z  is  obtained  by  taking  e  =  E  or  <  =  K  in  (8.6). 

The  tensile  and  shear  stresses  are  then  given  by  (8.4)  with 
these  expressions  for  £  and  Z  .  It  is  evident  that  they  are,  in 
general,  different  from  those  which  correspond  to  the  stepwise 
deformation  previously  considered.  Similar  disparities  will  evi¬ 
dently  be  found,  except  perhaps  in  some  exceptional  cases,  if  other 
forms  for  £  are  adopted. 

In  mathematical  terms,  the  type  of  behavior  predicted  by  the 
model  we  have  considered  arises  from  the  fact  that  the  intrinsic 
time  is  not  a  continuous  functional,  in  the  sense  of  the  supremum 
norm,  of  the  strain  history. 

We  note  that  if  a  material  does,  in  fact,  exhibit  the  type  of 
behavior  predicted  in  this  section,  it  will  be  extremely  difficult, 
and  perhaps  impossible,  to  subject  it  to  meaningful  tests  of  the 
type  considered. 
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9 .  Valanis's  second  theory  -  phenomenological  approach 

In  presenting  his  first  theory,  Valanis  claimed,  as  a  major 
advance,  that  the  endochronic  assumption  enables  us  to  construct 
a  continuum-mechanical  theory,  for  rate- independent  elastic- 
plastic  materials,  in  which  the  material  properties  are  described 
by  a  single  constitutive  equation  without  the  need  to  introduce  the 
concept  of  a  yield  surface.  These  constitutive  equations  also  have 
the  attractive  feature  that,  in  them,  the  stress  is  related  to 
total  strain. 

Both  of  these  features  are  lost  in  the  second  theory,  advanced 
by  Valanis  in  order  to  meet  some  of  the  criticisms  of  his  first 
theory.  This  second  theory  is  motivated  by  "thermodynamic"  argu¬ 
ments  of  a  type  similar  to,  and  open  to  the  same  criticisms  as, 
those  advanced  in  developing  the  first  theory. 

In  the  present  section  we  will  discuss  the  theory  in  more  con¬ 
ventional  phenomenological  terms.  In  [5]  Valanis  presents  his 
ideas  at  various  levels  of  generality,  but  here  we  will  limit  our 
discussion  to  that  form  of  the  theory  whose  implications  are  dis¬ 
cussed  in  [5]  at  greatest  length. 

The  strain  e  is  regarded  as  the  sum  of  an  elastic  strain 
Eg  and  a  plastic  strain  Ep  ,  thus: 

£  =  +  ep  .  (9.1) 

The  stress  a  is  related  to  the  elastic  strain  by  a  constitutive 
relation  of  the  form 

°  =  2yo£E  +  VtrEE^  ’ 


(9.2) 
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where  uQ  and  AQ  are  constants.  It  is  assumed  that  dilatational 
deformations  are  purely  elastic,  so  that 

tr  £p  *  Q  ,  tr  £  =  tr  £g  .  (9.3) 

From  (9.2)  and  (9. 3)2,  we  have 

tr  0  =  (2u0+  3AQ)tr  e  .  (9.4) 

From  (9.2),  the  deviatoric  stress  s  is  given  by 

s  -  a  -i(tr  a)  =  2yQeE  ’  (9.5) 

where  e^  denotes  the  deviatoric  elastic  strain.  From 

(9.1),  (9.3)  and  (9.5),  the  deviatoric  plastic  strain  ep  is  given 

by 

~P  =  ~P  =  ®  ‘  ®E  *  t  ‘  ?/2yQ  >  (9.6) 

where  e  denotes  the  deviatoric  total  strain. 

An  intrinsic  time  5  is  defined  in  terms  of  the  plastic 
strain  by 

d?  *  (trfdep)2}15  ,  (9.7) 

with  5  =  0  in  the  undeformed  state.  We  note  that  t,  is  the  arc 
length  measured  along  the  strain  path  in  the  9-dimensional  space 
defined  by  the  components  of  the  plastic  strain  in  a  rectangular 
cartesian  coordinate  system.'  [We  note  the  analogy  between  £  ,  so 
defined,  and  t  defined  by  (5.1).] 

For  a  rate -independent  material  the  plastic  strain  may  be 

regarded  as  a  function  of  5  .  It  is  assumed  in  [5]  that  the 
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deviatoric  stress,  when  ?  =  Z  ,  is  related  to  the  plastic  strain 
by  a  constitutive  equation  of  the  form 


s  =  s (Z) 


C=z 


<P(Z)  +  r  , 


(9.8) 


where 


r 


f(Z,OdepCO  , 


(9.9) 


is  a  positive  constant,  and  li'(Z)  is  a  positive  monotoni- 
cally  increasing  function  of  Z  .  From  (9.8)  and  (9.7)  it  follows 
that 


tr(s-r)2  =  [sYCo)KZ)]2  • 


(9.10) 


Equation  (9.10)  may  be  regarded  as  a  hypersphere  in  the  9- 
dimensional  space  formed  by  the  components  of  s  in  a  rectangular 
cartesian  coordinate  system.  The  center  of  this  hypersphere  is  at 
r  (regarded  as  a  vector  in  the  9-dimensional  space)  and  its 
radius  is  s^°-)^(Z)  .  Since  tr  s  =  0  ,  s  must  be  on  the  inter¬ 
section  of  the  hyperplane  tr  s  =  0  with  this  hypersphere.  This 
intersection  is  itself  a  hypersphere  in  the  8-dimensional  sub¬ 
space  of  the  9-dimensional  space  for  which  tr  s  =  0  .  This  hyper¬ 
sphere  is  called  the  yield  surface .  Since  tr  ep  =  0  ,  it  follows 
that  the  point  r  lies  in  this  subspace.  It  is,  of  course,  the 
center  of  the  hypersphere  in  the  8 -dimensional  space  and  the  radius 
of  this  hypersphere  is  s^°H(Z)  . 

We  note  from  (9.7)  that  if  dep  =  0  ,  then  d?  =  0  and, 
conversely,  if  d£  a  0  ,  then  dep  =  0  .  Thus,  if  at  any  instant 
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the  deformation  is  purely  elastic,  dep/d;^,,  indeterminate  ar*d 
equation  (9.8)  becomes  meaningless.  However,  from  (9.2),  the  stress 
increment  da  is  then  related  to  the  strain  increment  de(=deg)  by 

da  =  2uQde  +  (tr  de)<$  (9.11) 

and,  from  (9.5),  the  deviatoric  stress  increment  ds  is  related  to 
the  deviatoric  strain  increment  de(=de£)  by 

ds  =  2pQde  .  (9.12) 

If  deD/ds|  _  is  not  indeterminate,  it  follows  from  (9.8) 
that  the  plastic  strain  path,  at  arc  length  Z  ,  is  in  the  direc¬ 
tion  of  the  outward  normal  to  the  yield  surface  corresponding  to 
arc  length  Z  . 

We  see  that  the  theory  which  has  been  presented  is  of  a  type 
generally  similar  to  many  other  plasticity  theories  which  have  been 
formulated.  In  effect,  it  is  assumed  that: 

(i)  The  total  strain  may  be  regarded  as  the  sum  of  an  elastic 
strain  and  a  plastic  strain;  the  plastic  strain  is  isochoric. 

(ii)  It  is  assumed  that  the  moduli  associated  with  changes  of 
elastic  strain  are  constants,  independent  of  strain  history. 

(iii)  There  exists  a  spherical  yield  surface  in  deviatoric  stress 
space,  whose  outward  normal  is  parallel  to  the  incremental  plastic 
strain  vector, 

(iv)  The.  radius  of  the  yield  surface  and  the  position  of  its 
center  depend  on  the  history  of  the  plastic  strain  (i.e.  the  material 
modeled  exhibits  both  kinematic  and  isotropic  hardening) .  A  particular 
form  is  chosen  for  this  dependence,  which,  while  having  a  measure 
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of  generality,  is  far  from  being  of  the  most  general  form  that  can 
be  envisaged. 

(v)  If  the  deviatoric  stress  lies  at  a  point  on  the  yield  surface 
corresponding  to  some  specified  plastic  strain  history  and  is  then 
changed  to  a  point  lying  inside,  or  on,  this  yield  surface,  the 
corresponding  deformation  is  purely  elastic. 

It  is  evident  that  the  assumptions  (i) ,  (iii)  and  (iv)  meet 
the  objections  to  the  first  theory  of  Valanis  discussed  in  §  § 4 ,  5 
and  7.  Paralleling  the  discussion  in  §6,  we  shall  now  discuss  the 
extent  to  which  experiments  involving  superposed  simple  extensions 
and  shears  carried  out  in  discrete  time  intervals  can  be  used  to 
establish  the  validity  of  the  constitutive  equation  (9.8),  with 
(9.9)  . 

We  again  consider  a  somewhat  wider  class  of  constitutive 
equations  which  have  the  form  given  in  (9.8)  and  (9.9),  but  in 
which  ;  is  defined  in  a  different  manner.  Let  l p  denote  the 
length  of  the  plastic  strain  path  at  time  t  ,  thus  (cf.  (9.")) 

Ip  -  *p(t)  =  [*  (trtdcpU)]2}55  ,  (9.15) 

0 

and  let 

Lp  =  £p(T)  .  (9.14) 

Then,  Ip  can  be  used  to  parametrize  £p  . 

We  now  define  d t,  as  an  arbitrary  positive  isotropic  scalar 
function  of  dcp  ,  which  is  homogeneous  of  degree  unity  in  the 
latter.  Accordingly,  d£  is  a  positive  function  of  the  isotropic 
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invariants  of  dtp  ,  denoted  1^,  I 2,  I  and  defined  by  (cf.(5.6)) 
l±  =  tr  [dep  C-fp)  ]  =  0  ,  I2  -  tr  [dep  (£p)  ]  2 ,  I3  =  det  [dep  (£p)  ]  , 

(9.15) 

which  is  homogeneous  of  degree  unity  in  dep  .  We  take  5  =  0  , 
when  ip  =  0  .  We  define  Z  by 

Z  =  C (Lp)  .  (9.16) 

With  this  new  definition  of  £  ,  we  assume  a  constitutive 
equation  of  the  form  (9.8),  with  (9.9).  Then,  it  is  easily  shown, 
in  a  manner  analogous  to  that  employed  in  §6  in  discussing  the 
first  theory,  that  measurements  of  stress,  in  experiments  involving 
simple  extensions  and  shears  carried  out  in  disjoint  time  intervals, 
cannot  be  used  to  establish  the  validity  of  any  particular  form  of 
dependence  of  5  on  I2  and  1^  . 

Also,  it  can  be  shown,  paralleling  the  discussion  in  §8,  that, 
in  general,  for  a  constitutive  equation  of  the  form  given  in  (9.8) 
and  (9.9),  paths  in  plastic  strain  space  which  are  arbitrarily 
close  together,  may  have  very  different  values  for  c  correspond¬ 
ing  to  the  same  values  of  plastic  strain  and  accordingly  may  yield 
very  different  values  for  the  stress  at  the  same  value  of  the 
strain.  However,  this  fact  does  not  constitute  an  objection  to 
the  theory,  since  plastic  strain  histories  which  are  close  together, 
in  the  sense  of  the  supremum  norm,  may  be  associated  with  very 
different  histories  of  the  total  strain. 


Appendix  A.  "Thermodynamic"  justification  of  the  first  theory 
In  this  section,  we  outline  the  essential  features  of  the 
argument  of  Valanis  in  arriving  at  the  constitutive  equations  of 
his  first  theory  of  plasticity,  in  the  particular  case  when  the 
strains  are  small  and  the  deformations  are  carried  out  isother- 
mally.  The  argument  given  here  is  not  strictly  that  given  by  Vala¬ 
nis  [1].  His  theory  involves  rather  elaborate  thermodynamic  con¬ 
siderations  of  questionable  validity,  which  are,  in  any  case,  un¬ 
necessary  for  attaining  his  final  constitutive  equations.  We  have 
found  it  possible  to  replace  these  arguments  by  far  simpler  ones 
which  involve  only  some  of  the  assumptions  used,  either  explicitly 
or  implicitly,  by  Valanis  and  none  that  are  not  used  by  him. 

Valanis  defines,  for  a  rate  -  independent  material,  an  intrinsic 
time  scale  z  which  is  a  monotonically  non-decreasing  function  of  re 
time  t  and  is  a  functional^  of  the  strain-history  up  to  time  t  . 

He  considers  that  the  form  of  this  functional  dependence  may  depend 
on  the  material  considered.  For  an  isotropic  material  undergoing 
small  deformations,  he  expresses  z  as  a  monotonically  increasing 
function  of  another  (positive)  variable  5  ,  thus 

z  =  f  U)  .  (A .  1) 

S  is,  in  turn,  defined  by  an  expression  of  the  form 

(d s)2  =  k1  (tr  ds)2  +  k2tr(de)2  ,  (A. 2) 

where  k.  and  k2  are  positive  material  constants,  with  5=0 
initially.  We  note  that  the  expression  for  d£  is  the  most 


general  quadratic  form  in  da  ,  with  constant  coefficients,  which 
is  invariant  under  an  arbitrary  orthogonal  transformation,  Valanis 
then  assumes  that  the  stress  a  measured  at  time  t  is  an  iso¬ 
tropic  tensor  functional  of  the  strain  regarded  as  a  function  of  z  . 

These  assumptions  ensure  that  if  we  consider  the  strain  in  the 
material  to  execute  a  specified  path  in  nine-dimensional  strain 
space,  the  stress  corresponding  to  a  particular  point  on  this  path 
is  independent  of  the  rate  at  which  the  path  is  executed.  The  same 
objective  can  be  achieved  by  choosing  z  to  be  any  isotropic  scalar 
functional  of  the  strain  history  which  increases  monotonically  with 
time.  For  example,  the  strain  path  length  l  ,  defined  by  (5.2), 
provides  such  a  functional.  It  is  emphasized  that  the  particular 
choice  which  is  made  neither  increases  nor  decreases  the  generality 
of  the  theory  developed.  It  merely  changes  the  form  of  the 
functional  expression  for  a  appropriate  to  a  given  material. 

The  particular  definition  (A. 2)  of  £  adopted  by  Valanis  is 
justified  only  by  the  statement  [1]  "it  appears  logical  to  define 
?  by..."  .  The  particular  dependence  of  z  on  5  which  he 

assumes  is  avowedly  justified  only  by  alleged  agreement  of  the  pre¬ 
dictions  of  the  theory  with  experiment. 

From  this  point  onwards,  the  argument  of  Valanis  can  be  very 
much  simplified  without  losing  his  final  result.  We  shall  present 
this  simpler  argument  here  and  avoid  the  questionable  thermo¬ 
dynamic  argument  given  by  him  in  [1] . 

We  assume,  with  Valanis,  that  the  state  of  the  material  at 
any  instant  can  be  characterized  by  the  instantaneous  values  of  the 
strain  e  and  of  v  internal  variables  q ^ (a=l , . . . , v)  ,  which 
are  symmetric  second-order  tensors.  It  is  further  assumed,  still 
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r 


following  Valanis,  that  the  rate  of  change  of  q^a-*,  measured 

with  respect  to  z  ,  is  determined  by  the  instantaneous  values 

f  S'! 

of  e  and  of  q  (  3=1 ,  .  .  .  ,  v)  ,  thus 


dS 


(cO 


Iz 


£  (CC) 


)  . 


(A.  3) 


fa'I 

This  is  the  evolution  equation  for  q  1  .  (We  will  see  later 

that  this  assumption  is  unacceptable,  but  for  the  moment  will 

pretend  that  this  is  not  the  case.) 

(a) 

The  function  J  must  be  an  isotropic  symmetric  tensor 

function  of  e  and  q^®^,  since  the  material  is  isotropic.  We 

shall  suppose  that  it  depends  sufficiently  smoothly  on  e  and 

q*-®-*  and  that  these  are  sufficiently  small,  so  that  we  can  neglect 

terms  of  higher  degree  than  the  first  in  them.  We  also  make  the 

fa'I 

assumption  that  s  =  qv  }  *  0  when  z  =  0  .  We  accordingly  write 


dq 


(a) 


dz 


=  A  s  + 

a~ 


v 

I 

S=1 


Aas3(8)  +  ^Batr~  + 


v 

8=1  p 


tr  q(S))5 


(A.  4) 


where  the  A's  and  B's  are  constants  and  o  is  the  Kronecker 
delta.  From  (A. 4)  we  readily  obtain 


d(tr  v 

—  =  (A  +5B  )tre  +  7  (A  +  3B  -)tr  q 

v  a  af;  v  a8  aS'  ~ 

P“  i- 


■ar 


(3) 


,  (a) 

d9o 


(A.  5) 


3z  =  Aa?D  +  Aa89D 


CS) 


where  the  deviatoric  strain  and  deviatoric  internal  variables  are 
given  by 
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:d  =  £  -»(tre)5  ’  9d^  =  9^  "J^tr  • 


(A.  6) 


r  ct  i 

Equations  (A.5)?  can  be  solved  for  to  yield  a 

solution  of  the  form 


5°“’  *  eL  *»-> 


V  _  [Z  -pft(z-z') 

e  3  5D(2,)dz 


(A. 7) 


0 


where  the  A's  are  constants  determined  by  the  equations 


A3v> ^Aa8+pp5afP  ~  0  ’ 

3-1 


^  Aa$  Aa  ’ 


6=1 


(a ,y=l , . . . ,v) 


(A.  8) 


in  which  p  (y=l,...,v)  are  the  solutions  for  p  of  the 
equation 


det|p6a8  +  AaS |  =  0  ’ 


(A.  9) 


where  6  „  is  the  v-dimensional  Kronecker  delta. 
aB 

In  a  similar  manner  equations  (A. 5)^  can  be  solved  for 
tr  q*-a-*  to  yield 


tr  q 


(*)  = 


l  B 


6  =  1 


a6 


-\Az-z  ) 

f  *  tre(z’)d2' 


(A. 10) 


where  the  B's  are  constants  given  by 


l_  B3u(Aa3+3Ba8  +  Xu6a6^  =  0  » 

f-1  =  1 

(a,u=l,...,v)  (A. 11) 


l  B. 


»  A, 


3B. 


8=1 
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in  which  X  (u=l . v)  are  the  solutions  for  X  of  the 

equation 


det I X6 


+  3B 


(A. 12) 


In  deriving  (A. 7)  and  (A. 10),  we  make  the  assumption  that 
=  0  when  2  =  0. 

With  the  assumption  that  s(z)  =  0  when  2=0,  we  can 
rewrite  (A. 7)  and  (A. 10)  as 


L,  1  f2  '’bC*-*') 


den(z')l, 


(A. 13) 


tr  3*a)  “  l  Ba ( y~  tre(z)  -p- 

3=1  a“  vAS  ~  Ag 


^  -Xft (z -z ' ) 


dltrefz  ')] 


We  now  assume  that  the  Cauchy  stress  a  at  "time"  z  is  an 
isotropic  linear  function  of  e(z)  and  q^a)  (a=l,...,v)  .  The 
deviatoric  stress  aQ  is  defined  by 

2d  =  H  "  I  •  (A. 14) 


Then,  we  may  express  aQ  and  tra  in  the  forms 

aD  =  G2d  +  G32d^  5 

tra  =  H(trs)  +  £  H.(tr  q^)  , 
S=l  d  ~ 


(A.  15) 


where  the  G’s  and  H’s  are  constants.  Using  (A. 13),  we  obtain 
from  (A. 15),  expressions  for  aQ  and  tra  of  the  forms 
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?D 


p(2-Z  ' )d£p(Z  '  )  , 


tra 


k(z-z ’)d[trs]  , 

0 


(A. 16) 


where  the  functions  y(z-z') 


and  k(z-z')  have  the  forms 


y (z-z  ' ) 
< (z - z  '  ) 


y  -pfi(z-z') 

V 


<Q  +  l  <Be 
0  8=1  * 


-XQ(z-z') 


(A. 17) 


and  Hg,  <g  (8=0,..., v)  are  constants. 

It  is  evident  that  the  passage  from  equations  (A. 13)  and 
(A. 15)  to  the  constitutive  equations  (A. 16)  is  valid,  whatever  the 
physical  interpretation  of  z  ,  provided  that  it  is  an  isotropic 
variable  in  terms  of  which  the  histories  of  the  strain  and  of  the 
internal  variables  can  be  parametrized.  The  choice  of  z  made 
by  Valanis  is 

z  =  £  (?)  ,  £(?)  =  |  £n(l  +  SO  ,  (A. 18) 


where  3  is  a  positive  constant  and  c  is  defined  by  (A. 2). 

If  we  accept  a  relation  of  the  form  (A.18)^  and  use  (A. 3), 
we  obtain 

dq 

-jt—  =  £’  U)f(a)  U>q(s))  •  (A. 19) 

This  implies  that  the  infinitesimal  changes  in  the  internal 
variables  due  to  an  infinitesimal  change  in  the  strain  depends  not 
only  on  the  instantaneous  values  of  the  internal  variables  and  the 
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strain,  but  also,  through  ;  ,  on  the  whole  past  history  of  the 
strain.  We  remark  that  usually,  when  internal  variables  are  intro¬ 
duced,  along  with  the  current  strain,  as  independent  variables  in 
a  constitutive  equation,  it  is  in  order  to  provide  a  full  des- 
cription  of  the  current  state  in  terms  of  the  current  values  of  the 
independent  variables  of  the  theory.  If  the  internal  variables  and 
the  strain  provide  a  complete  description  of  the  state,  then  the 
infinitesimal  change  in  the  internal  variables  due  to  a  specified 
infinitesimal  change  in  the  strain  should  depend  only  on  their 
current  values. 

Another  peculiarity  of  the  relation  (A. 19)  and,  indeed,  of 
(A. 3)  is  that  different  infinitesimal  changes  of  strain  which  result, 
from  (A. 2),  in  the  same  values  of  dg  ,  lead  to  the  same  infini¬ 
tesimal  changes  in  the  internal  variables  .  This  peculiarity 

can  be  avoided  by  including  the  "rate"  of  change  of  strain,  de/dz  , 
as  an  independent  variable  in  (A. 3).  Then,  the  relations  (A. 5) 
are  replaced  by 


d(tr  q^a)) 
dz 


d(tre) 

(A  +3B  ) tre  +  (C  +3D  )  — ,-C- 
a  or  v  a  or  dz 


dz 


■lx  S'” 


dcD  V 


Aa?D  +  Ca  JT+  l  Aas9D6)  ’ 

3=1 


(A. 20) 


where  the  A's,  B's,  C's  and  D's  are  constants. 

From  (A. 16),  proceeding  as  before,  we  again  obtain  expressions 
for  and  tra  of  the  forms  (A. 16)  where  y(z-z')  and  k(z-z') 


still  have  the  forms  (A. 17),  and  u,  and  <g  (3=1 , . . . ,v)  are 
(different)  constants,  o,  and  \a  (3=1,..., v)  are  still  the 
solutions  of  equations  (A. 9)  and  (A. 12).  We  note,  however,  that 
in  obtaining  (A. 17)  from  (A. 20)  we  must  either  assume  that 
f^(  )  in  (A. 3)  depends  linearly  on  de/dz  ,  or  that  it  is  a 
sufficiently  smooth  function  of  de/dz  and  that  de/dz  is 
sufficiently  small.  The  assumption  that  e  and  J  are  small 
does  not,  of  course,  guarantee  that  de/dz  is  small.  Indeed, 
de/dz  will  not,  in  general,  be  small.  For  example,  if  the  defor¬ 
mation  considered  is  a  simple  shear,  then  the  shear  component  of 
-  h 

de/dz  is  C )  ,  while  the  remaining  components  are  zero. 

Even  when  modified  by  the  inclusion  of '  de/dz  as  an 
indeoendent  variable,  the  evolution  equation  (A. 19)  still  implies 
that  infinitesimal  changes  in  the  internal  variables  resulting 
from  a  specified  infinitesimal  change  or  strain  is  indepencent 
of  whether  the  strain  is  purely  elastic  or  plastic.  It  is 
difficult  to  see  what  phys ical  identification  of  the  internal 
variables  could  lead  to  such  a  result. 
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